Particle acceleration in supernova-remnant shocks 



S.P. Reynolds 



Abstract It has been known for over 50 years that the 
radio emission from sheU supernova remnants (SNRs) 
indicates the presence of electrons with energies in the 
GeV range emitting synchrotron radiation. The dis- 
covery of nonthermal X-ray emission from supernova 
remnants is now 30 years old, and its interpretation 
as the extension of the radio synchrotron spectrum re- 
quires electrons with energies of up to 100 TeV. SNRs 
are now detected at GeV and TeV photon energies as 
well. Strong suggestions of the presence of energetic 
ions exist, but conclusive evidence remains elusive. Sev- 
eral arguments suggest that magnetic fields in SNRs are 
amplified by orders of magnitude from their values in 
the ambient interstellar medium. Supernova remnants 
are thus an excellent laboratory in which to study pro- 
cesses taking place in very high Mach-number shocks. 
I review the observations of high-energy emission from 
SNRs, and the theoretical framework in which those 
observations are interpreted. 



1 Inferences from radio emission 

Remnants of historical supernovae had been known 
since Lundmark's identification of the Crab Nebula 
with SN 1054 AD (jHubbldl 19281) . The first radio source 
to be identified as a previously unknown supernova rem- 
nant (SNR) was Cassiopeia A, proposed by Shklovskii 
in 1953, who also suggested that synchrotron radia- 
tion was the mec hanism for producing radio emission 
(jShklovskiil Il953l ). based on the observed power-law 
spectrum Si, oc with a ^ 0.8 for Cas A. Minkowski 
in 1957 confirmed the ide ntification from optical obser- 
vations (| Minkowskil 1 1 9 5 7l) . Thus from the early days of 
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radio astronomy, it was recognized that energetic par- 
ticles were present in SNRs. Basic synchrotron physics 
tells us that emission at frequency v is primarily by elec- 
trons with energy E = 15(i^(GHz)/B(^G))i/2 QeV, so 
radio emission at frequencies of a few hundred MHz, 
typical at the time, immediately implied the presence 
of electrons with Lorentz factors of 10"^ — 10'*. 

Radio remains the spectral region in which SNRs 
are most consistently identified. Green's famous cata- 
logue of Galactic SNRs (Green 2009; available online at 
http: //www.mrao. cam. ac.uk/surveys/snrs/l lists 274 
objects, essentially all with known radio properties. 
The mean spectral index is about 0.5, but with a sig- 
nificant spread of order 0.2 (Figure [I]). This spread is a 
significant problem for theories of particle acceleration 
described below. Also significant is a tendency for the 
historical remnants (less than 2000 years old) to have 
steeper indices (a > 0.6), a trend continued by radio su- 
pernovae, not represe nted here, which ca n have indices 
as steep as 0.9 - 1.0 (jWeiler et al.l 120091 ). A few rem- 
nants such as Cas A have very well-sampled radio spec- 
tra (see references in Green 2009); Cas A's spectrum is 
well described by a single power-law with spectral index 
of 0.77 between 100 MHz and 100 GHz. However, most 
remnants are represented by only a few data points 
with substantial error bars. Other historical remnants 
have spectra with suggest ions of concave-up curvature 
(I Reynolds fc Ellisonlll992l ). naturally explained by effi- 
cient shock acceleration (see below). 

Magnetic field strengths are difficult to measure in 
SNRs; since the intensity of synchrotron radiation from 
a power-law distribution of electrons N{E) — KE^^ 
electrons cm^'^ erg^^ scales as KB^^", radio syn- 
chrotron fluxes only give roughly the product of the 
energy densities in magnetic field and electrons. X- 
ray evidence described below indicates that magnetic 
fields are substantially amplified over typical interstel- 
lar values of a few microGauss, but the only avenue 
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for estimating energetics from radio data is the use of 
minimum-energy equipartition arguments. These ar- 
guments indicate that the minimum energy in elec- 
trons and magnetic field in typical SNRs is far below 
the ^ 10^^ erg explosion energies; SNRs are not ef- 
ficient synchrotron radiators. The equipartition mag- 
netic fields so derived tend to be low, but as there is 
no strong physical argument that equipartition should 
hold, or even among which particles (should one in- 
clude ions?), the equipartition magnetic field strength 
is really just a proxy for mean surface brightness. 

Magnetic-field orientations can be usefully derived 
from radio polarization directions. A uniform syn- 
chrotron source with spectral index a — 0.5 has a 
polarized fraction of about 70%, but very few rem- 
nants show values above 40%. Young SNRs have much 
lower values, typically of order 10% - 15% (see ref- 
erences in Reynolds & Gilmore 1993), implying that 
their magnetic fields are primarily disordered. The or- 
dered components, however, tend to be radial. In older 
remnants, magnetic-field orientations are typically con- 
fused, but it is fairly common to see a tangential ori- 
entation, which one would expect if a high-compression 
radiative shock compresses upstream magnetic field, in- 
creasing the tangential component by a factor of the 
compression ratio. 

Even though radio observations of SNRs do not re- 
quire a large fraction of the SN energy, it is possible to 
argue that young SNRs require acceleration of new elec- 
trons - simply borrowing and compressing relativistic 
electrons from the cosmic-ray population in the ISM 
is inadequate both because of high observed surface 
brightnesses and spectra much different from tho se of 
low-energy cosmic-ray electrons (jRevnoldsl l2008al) . 

Radio observations of SNRs have left several unex- 
plained puzzles, most several decades old now. What is 
responsible for spectral indices flatter than 0.5? (While 
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a few pulsar-wind nebulae contaminate the low-a end 
of the distribution of Figure 1, most of the remnants 
with a < 0.5 are shells.) In a few cases, contamination 
with flat-spectrum thermal emission may be responsi- 
ble. Steeper spectra than 0.5 can be obtained with very 
low Mach-number shocks, but this requires < 10 
and is unlikely to be the case for as many remnants 
as Figure [1] requires. Finally, the question of the ra- 
dial orientation of the ordered component of magnetic 
field in young remnants remains unexplained, though 
the operation of fluid instabilities at the contact inter- 
face between shocked ISM and shocked ejecta is often 
invoked (e.g., Jun, Jones, & Norman 1996). 



2 X-ray synchrotron emission 

The featureless X- ray spectrum of th e remnant o f SN 
1006 above 1 keV |Becker et al.lll98d ITooiIIiQSOI ) was 
first attributed to sync hrotron radiation frorn shock- 
accelerated electrons by iRevnolds fc Chevalieii (jl98ll) . 
This explanation required the presence of electrons with 
TeV energies. However, detection o f oxygen lines from 



the spatially integrated spectrum (jGalas et al.l 119821 ) 



indicated the operation of thermal processes. Only 
with the ASCA observation of thermal line emission 
from the interior and a featureless continuum from the 
bright limbs (Koyama et al. 1995) was it clear that a 
synchrotron component was called for, though it repre- 
sented the rolling off of the spectrum due to some pro- 
cess limiting the maximum energy of electro ns, rather 
than to a straight extrapolation from radio (jRevno Idi 
I199d). 

How do wc know the featureless continuum emis- 
sion is synchrotron radiation? Nothing else works. A 
nonthermal power-law electron distribution with en- 
ergies of order tens of keV will produce a power-law 
bremsstrahlung spectrum with photon energies smaller 
by a factor of a few. However, those electrons would 
be just as efficient at exciting atomic lines as electrons 
drawn from a Maxwellian distribution. A tiny corner of 
parameter space might be available for a plasma either 
not yet ionized to the stages required for X-ray lines, o r 
on the other hand totally stripp ed (e.g 



LaminJl998l) 



Fig. 1 Spectral-index distribution of SNRs in Green's cat- 
alogue with adequately measured radio spectra, in bins of 
0.1 in Q. 



However, the careful study by iLamina (|1998[ ) shows 
that in detail this cannot work for SN 1006. Finally, 
inverse-Compton upscattering of any known source of 
photons would produce a spectrum with the same slope 
as that of the synchrotron emission produced by those 
electrons - that is, the radio slope (a ^ 0.6), far to o 
flat for the observations {ax ^ 2.3; Long et al. 20031) . 
Furthermore, the synchrotron hypothesis makes a pre- 
diction: spectral steepening to higher X-ray energies. 



3 



For SN 1006, this ha s been confirmed by INTEGRAL 



( Kalemci etallbood) . 



There are now four known Galactic remnants whose 
soft X-ray spectrum is dominated by synchrotron X- 
rays: i n addition to SN 1006, Gl. 9+0.3 dRevnolds et al. 



2008bl): G347.3- 0.5 (also known as GX J17 13. 7-3946 



after each return, but a decreasing number remains for 
further cycles, producing a power-law distribution (see 
Bell 1978 for a kinetic-theory probabilistic argument). 
A well-known result of this process, for energetically 
unimportant test particles, is a power-law distribution 

^ with index dependent only on the shock compression 

Slane et al.lll999l): and G266.2-1.2. "Vela Jr." (|Aschenbach| -atio: N{E) cx E-' with s = {r + 2)/{r - 1) where 

r = P2/P1 is the shock compression ratio, equal to 
4 for strong shocks. This result applies to extreme- 
relativistic particles {E = pc, p the momentum and 
c the speed of light). The prediction is thus s = 2 
which, for electrons, implies a synchrotron spectral in- 
dex a = (s — l)/2 = 0.5, in tolerable agreement with 
spectra observed from Galactic SNRs. 

The maximum energy to which particles can be ac- 
celerated depends on the limiting mechanism. We as- 
sume a diffusion coefficient k which scales with particle 
energy. This results if the particle mean free path is 
a multiple rj of its gyroradius: Amfp = Wg = rjE/eB 
for ultrarelativistic particles, since = jmc^ /eB (cgs 
units). Then the "Bohm limit" in which the mean free 
path is a gyroradius is 77 = 1. For weak turbulence, one 
expects 1] > 1, though this may not be a hard physical 
limit. For a remnant of age t with shock speed Ush, with 
surroundings containing MHD scattering waves only up 
to a wavelength Amax, the maximum energies scale as 



1998 t ISlane et al.ll200ll) . (See Reynolds 2008a for a 
more detailed review of SNR X-ray and 7-ray emis- 
sion.) Figure [2] shows the Chandra images of SN 1006 
and G1.9-f0.3. For both objects, thermal lines have 
been detected from fainter regions of the remnant. The 
low-energy (central) emission in Fig. [2] is primarily 
ejecta emission dominated by oxygen. Figure |4] shows 
the spectrum from the interior of Gl.9+0.3, with clear 
emission lines of Si, S, and Ar. 

Synchrotron X-ray emission contributes to the spec- 
trum of several more Galactic remnants (reviewed in 
Reynolds 2008a) . In historical (or quasi- historical) rem- 
nants RCW 86 (SN 185?), Tycho (SN 1572), Kepler 
(SN 1604), and Gas A (SN - 1680), "thin rims" of fea- 
tureless X-ray emission lie at the edges of the remnants, 
and are presumed to indicate the outer blast wave. Fur- 
thermore, all these remnants along with SN 1006 have 
been reported to have hard X-ray continua in integrated 
spectral (non-imagi ng) observations with the PGA in- 
strument on RXTE (|Allen. Gotthelf. fc Petrll999l ). It 
appears that synchrotron X-ray emission from the blast 
wave is a common feature in young SNRs, less than a 
few thousand years old. The implications are strik- 
ing: electrons are present in these objects with ener- 
gies E = 72{hv/l keV)i/2(s/io pG)-^/"^ TeV. Now the 
electron distributions are not straight power-laws all the 
way from radio-emitting energies (ca. 10 GcV); rather, 
in all known c ases, the observed X-rays fal l below that 
extrapolation ( Revnolds fc Keohanel (l999l ). indicating 
that some limitation on electron energies is taking place 
below 100 TeV (perhaps far below, for those older SNRs 
with no evidence for synchrotron X-ray emission). As 
will be summarized below, shock acceleration may be 
limited by the finite age (or size) of the remnant, par- 
ticle escape above some energy due to absence of scat- 
tering waves upstream, or (affecting electrons only) ra- 
diative losses. 



3 Diffusive shock acceleration 

In stand ard diffusive shock acceler ation (amply re- 
viewed in iBlandford fc Eichleil E987I ) . particles scatter 
from magnetic inhomogeneities borne by converging flu- 
ids on either side of a shock wave, gaining energy with 
each reflection. Most particles disappear downstream 
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In all cases, for Ush > 1000 km s~^ and ages above a 
few hundred years, maximum energies of 10 - 100 TeV 
are easily obtainable. The accelerated-particle spectra 
should show exponential cutoffs with these flducial en- 
ergies. For synchrotron emission, the dependence of 
the peak frequency iy„i on electron energy of Vm oc E^ 
means that the synchrotron spectrum will then drop 
roughly as , that is, considerably slower 

than exponential, and hardly differing from a power-law 
in the bandpass of most X-ray observatories (^ 0.3— 10 
keV). 

The diffusion coefflcient may be anisotropic; in par- 
ticular, diffusion along and across magnetic-field lines 
is likely to take place at different rates, with effects on 
the acceleration time r to some energy. If the shock 
velocity makes an angle ^Bn with the mean upstream 
magnetic field, we can parameterize this effect with 
Rj{0Bn,Vir) = r (6*311 )/t (6*311 = 0). We scale to typ- 
ical values for young SNRs: ug.5 = Ush/3000 km s~^; 
ts = t/1000 yr; Bio = S/10 ^G; and A17 = A„iax/10" 
cm. The frequencies at which electrons with energy 
Emax emit their peak power for each cutoff mechanism 
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are then 
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Of course, in a given object, the lowest value of i^max 
will be the operative value. Thus if one can determine 
the mechanism causing the spectral cutoff, its value 
constrains considerably more physical parameters than 
the simple observation of a radio power-law spectrum. 



models for these f our ca ses have been constructed (e.g., 
Berezhko fc Volkl l2006f ). The TeV emission in these 
models can be due either to IC from the CMB, or to 
tt" decay. Both classes of model have difficulty. The 
former require low filling factors of magnetic field, and 
imply inefficient shock acceleration, while the latter suf- 
fer from severe limits on thermal gas from X-ray obser- 
vations, implying insufficien t targets for the relativistic 
protons feU ison et al. 2010l ) . For complex objects such 
as G347.3-0.5 and Vela Jr., simple one-zone models may 
be inadequate, though there is as yet no clear path to 
consistent models. 



4 Radiation from GeV to TeV photon energies 

Populations of relativistic ions and electrons can pro- 
duce observable continuum radiation through four 
mechanisms, one hadronic and three leptonic (reviewed 
in Reynolds 2008a). The hadronic mechanism is the 
inelastic scattering of cosmic-ray protons on thermal 
nuclei, producing pions. The charged pions decay to 
electrons and positrons, making a (probably) negli- 
gible contribution to the relativistic lepton pool, but 
the TT^'s decay to gamma rays of comparable energy 
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70 MeV. The spectrum of 



emitted photons should be that of the ions that pro- 
duce them. The three leptonic processes are syn- 
chrotron emission, described above, as well as nonther- 
mal bremsstrahlung, with the same spectrum as that 
of the nonthermal electrons, and inverse-Compton (IC) 
upscattering of photons from any significant ambient 
radiation. In practice, this is likely to be primarily 
the cosmic microwave background (CMB), though in 
some cases, IC from UV-optical-IR photons may be 
competitive. The spectrum will be the same as that of 
the synchrotron emission from whatever population of 
particles is responsible, that is, considerably harder in 
the keV - TeV range than that of the other processes. 
While the synchrotron process is clearly operating, it 
is not clear which of the other processes might be re- 
sponsible for emission from any particular object. The 
best evidence for ion acceleration is the spectral fea- 
ture resulting from the minimum photon energy from 
a created pion nea rly at rest, about 7 MeV. Detailed 



calculations fe.g.. iBaring et al.l 119991 ) show that this 



feature may not be highly distinct in a real object. 

Several shell (i.e., not containing a pulsar) SNRs 
have been detected in TeV gamma rays, using air- 
Cerenkov detectors such as the High-Energy Stereo- 
scopic System in Namibia. These include G347.3-0.5, 
Vela Jr. RCW 86, and SN 1006 (see references in 
Revno iS l2008al) . The spectra are steep, with photon 
indices T - 2 (F^ cx {hi^yr); for G347.3-0.5, the spec- 
trum is observed to steepen above 1 TeV. Elaborate 



5 Magnetic-field amplification 

An important result first suggested many years ago, but 
only recently put on a firmer observational foundation, 
is the increase in magnetic-field strength in SNRs over 
a simple factor of the compression ratio expected in 
highly ionized gases. The possibility of much stronger 
shock amplification of magnetic fields was first pro- 
posed for helio spheric shocks by Chevalier ( 1977), while 
the models of Revnolds fc Chevalier ( 198ll ) demanded 
substantial amplification. Early gamma-ray upper lim- 
its from Cas A bounded the electron population from 
above (from the inferred absence of bremsstrahlung), 
bounding the magnetic-field strength from below based 
on obs erved radio synchrotron fluxes. ICowsik fc Sarkai 
( 1980l) used this argument to deduce a minimum mag- 
netic field strength of about 1 niG. The slight concave 
curvature observed in radio spectra of SNRs was ex- 
plained bv lRevnolds fc EllisonI (119921 ) as due to efficient 
shock acceleration modifying the shock structure, but 
values of magnetic field were also implied of 100 /iG 
and more - far higher than a few times the typical in- 
terstellar magnetic field of 3 - 5 /iG. 

High-resolution X-ray images from Chandra show 
that the "thin rims" present in most historical shell 
SNRs are in fact very thin ~ so thin that an unusual 
depletion (beyond simple post-shock expansion) of ei- 
ther relativistic electrons or magnetic field is required to 
explain the sudden disappearance of synchrotron emis 



sivity downstream (jBamba et al.ll2003l : IVink fc Laming 
20031 ). If synchrotron losses a re depleting the electrons. 



we infer ( Parizot et al. 20061 ) 



B > 200w8/^(w/0.01 pc)-^/VG 



2/3, 



(7) 



where the shock speed us = m/1000 km s~-^ and the 
filament width is w. However, if the amplified mag- 
netic field is in waves of some kind, these waves might 
damp, providing an alternative mechanism for the rims 
( Pohl et al.l l2005f ) . Figure [3] shows a typical rim - but 
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it is present in radio as well, where synchrotron losses 
cannot possibly be a factor. It is possible that both 
mechanisms are required to explain some rims. 

An additional argument for the amplification of 
magnetic field comes from the variations seen in pe- 
riods of a few years i n X-ray feat ures in G347.3-0.5 
(Uchivama et al.l 120071) and Gas A ( Patnaude fc Fesen 



20071) . Both brightening and fading are seen. If the rel- 



evant timescales are those of particle acceleration and 
synchrotron losses, in either case one obtains values of 
the magnetic field of order 1 mG. However, if the post- 
shock magnetic field is highly turbulent, one expects 
such "twinkling" due t o fluctuations even if the mean 
magnetic field is lower ( Bvkov et al. 20081) . 



6 The youngest Galactic remnant Gl.9+0.3 



The discovery of a remnant only about 100 yr old, from 
expansion between radio observations made in 1985 and 



Fig. 2 Chandra imag es of (top) SN 1006 ( NASA/CfA) 
and (bottom) Gl.9-1-0.3 (|Revnolds et al.ll2008bl ) bet ween 1.5 
and 7 keV, smoothed with platelets as described in IWilletd 
(l2007l ). 




X-ray Chandra observations in 2007 (jRevnolds et al 



2008bl) . opens a new window on the very early develop 



ment of a SNR, and of the physics of unprecedentedly 
fast SNR shocks. At an assumed distance of the Galac- 
tic Genter, 8.5 kpc (the very high absorbing column, 
Nh ~ 6 X 10^^ cm~^, means the distance cannot be 
much less), the mean expansion speed is about 14,000 
km s~^. The integrated synchrotron spectrum has a 
rollofF frequency of hiy = 2.2 keV, one of the highest ever 
measured. A recent longer observation has revealed the 
presence of spectral lines from the radio-bright N rim 
and interior (Figure S]) . 

The figure shows clear lines of helium-like states of 
Si, S, and Ar from the central region. Ga, while visible 
from the N rim, is absent in the center ~ but a line at 4.1 
keV, which we identify as due to ''^Sc, is clearly present. 
^''Sc is produced by electron capture in ■*'*Ti, a radioac- 
tive element with mean life of 89 yr. The inner-shell 
vacancy is filled by the emission of a 4.1 keV photon. 
This is the first firm detection of this element in a SNR. 





Fig. 3 Top: Chandra image of Tycho's SNR (NASA /CfA). 
Below: VLA image at 1.4 GHz (|Reynoso et al.lll997^ . Note 
the thin rim in both bands at the NE edge (upper left). 



Fig. 4 Spectrum from the central region of Gl.9+0.3 
(Borkowski et al., in preparation). 
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The hne strength we measure imphes an initial mass of 
^''Ti of (1 - 3) X 10"'^ Af© for an age of 100 - 140 
yr, within the predicted range for either core-collapse 
or Type la supernovae, though earlier spherically sym- 
metric Type la models predict lower 44^j masses. Since 
the ''^Sc need not be ionized to produce the 4.1 keV X- 
rays, we are sensitive to both shocked and unshocked 
'^'^Sc. A longer observation can allow the determination 
of the spatial distribution and kinematics of ^^Ti, of 
considerable interest to SN modelers. 



7 Status report 

1. X-ray observations show that SNR shocks routinely 
exhibit synchrotron radiation from electrons up to 
100 TeV in energy. 

2. Standard diffusive shock acceleration can easily ac- 
count for these energies for a range of conditions. 
Little, however, can be predicted. Quantities of in- 
terest include efficiencies, obliquity-dependence (i.e., 
^Bn-dependence), or e/p ratio. 

3. Emission at TeV energies seen from several young 
SNRs may be leptonic (IC/CMB) or hadronic (tt"- 
decay). Both types of models have serious problems. 
In particular, while many indirect lines of evidence 
suggest that these shocks efficiently accelerate ions, 
direct observational confirmation is still lacking. 

4. Both integrated fluxes and small-scale structure in 
X-rays indicate that magnetic fields in SNRs are 
higher by 10 - 100 than expected from simple com- 
pression of ISM fields. This may b e due to a cosmic- 
ray driven instability ( Belli 2004 ). This instability 
requires high-efficiency ion acceleration. 

5. Future needs: Observations at GeV energies (and, 
we hope, down to 100 MeV or below) with the Fermi 
Gamma-Ray Space Telescope will be important in 
constraining models. Theoretical work to predict 
how particles are injected into the shock-acceleration 
process as well as the items enumerated above is 
essential. Radio observations provide information 
which has still not been used by modelers to maxi- 
mum effect. Finally, Gl. 9-1-0. 3 is evolving so quickly 
that its monitoring with time, and the accumulation 
of better statistics for spectroscopy, may add sub- 
stantially to our understanding of shock acceleration 
in young SNRs. 
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